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ABSTRACT 

The inferred black hole in the Galactic center spans the largest angle on the sky among all 
known black holes. Forthcoming observational programs plan to localize or potentially resolve 
the image of Sgr A* to an exquisite precision, comparable to the scale of the black hole horizon. 
Here we show that the location of the image centroid of Sgr A* should depend on observing fre- 
quency because of relativistic and radiative transfer effects. The same effects introduce a generic 
dependence of the source polarization on frequency. Future detection of the predicted centroid 
shift and the polarization dependence on frequency can be used to determine the unknown black 
hole spin and verify the validity of General Relativity. 

Subject headings: Galaxy: center — submillimeter — infrared: general — black hole physics — accretion, ac- 
cretion disks — gravitational lensing 



1. Introduction 

The individual orbits of fast-moving stars near 
the dynamical center of the Milky Way galaxy 
indicate the existence of a black hole of mass 
~ 4 X 10^ solar masses at the position of the ra- 
dio/infrared source Sgr A* (Schodel et al. 2003; 
Ghez et al. 2005). The angular scale occupied 
by the horizon of this black hole, 5-10 micro- 
arcseconds (/^as), is a factor of ~ 2 larger than the 
nuclear black hole in the distant galaxy M87 and 
orders of magnitude larger than all other known 
black holes. Recent observations revealed infrared 
emission from the radio source Sgr A* which flares 
on a timescale of tens of minutes, comparable to 
the horizon crossing time of accreting gas (Gen- 
zel et al. 2003). This data has motivated an ob- 
servational program to monitor shifts of tens of 
jxas in the centroid of the infrared image of Sgr 
A* due to an orbiting hot spot in the accretion 
flow during flares (Broderick & Loeb 2005), using 
the newly constructed PRIMA instrument on Very 
Large Telescope (VLT) (Paumard et al. 2005). 
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In parallel, there are plans to image Sgr A* at 
~ 20/xas resolution with a Very Large Baseline 
Array (VLBA) of observatories at sub-millimeter 
wavelengths (Miyoshi et al. 2004), at which the 
accreting gas is expected to be optically-thin and 
the black hole shadow may be recovered (Falcke 
et al. 2000). 

The faint emission from Sgr A* is commonly 
interpreted as a consequence of a radiatively inef- 
ficient accretion flow around the black hole (Yuan 
et al. 2003). The observed break in the sub- 
millimeter spectrum is thought to signal the tran- 
sition of the flow to the optically thin regime at 
wavelengths below the break. The flux below the 
break wavelength is expected to be dominated by 
emission from the vicinity of the black hole hori- 
zon, where general relativistic effects become sig- 
niflcant. These effects result from the warping 
of spacetime by the black hole (e.g., gravitational 
redshift, rotation of polarization via parallel prop- 
agation) and the rapid motion of the disk near 
the innermost stable circular orbit (ISCO; e.g., 
Doppler shift, beaming, Lorentz boosted polariza- 
tion vectors, etc.). Although attempts have been 
made to incorporate these effects into the deter- 
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mination of the unpolarizcd spectrum in the opti- 
cally thin limit (Kurpiewski & Jaroszynski 2000), 
earlier studies have not included a complete treat- 
ment of the polarization spectrum^ and are not ap- 
plicable to the sub-millimeter wavelength regime 
where the opacity is non-negligible. 

Here we produce images for three variations 
on the model of Yuan et al. (2003), listed in Ta- 
ble 1, each adjusted to roughly fit the observed 
radio and infrared spectrum. The primary dis- 
tinction between the different cases in Table 1 
is the black hole spin, ranging from non-rotating 
to maximally rotating (a = 0.998). This was 
done using the methods employed by Broderick & 
Loeb (2005) and detailed in Broderick & Bland- 
ford (2003, 2004), and we direct the reader to these 
papers for more information. From the calculated 
images we compute the unpolarized and polarized 
spectra as well as the location of the image cen- 
troid as functions of frequency. 

2. Accretion Models 

We model the structure of the accretion flow 
based on the results of Yuan et al. (2003), who 
showed that the vertically-averaged electron den- 
sity and temperature could have a nearly power- 
law dependence on radius. Here we write the den- 
sity, rig) and temperature, Te, of the thermal elec- 
trons, as well as the density, rinth, of non- thermal 
electrons, as having the spatial distribution 

„ / r \ -0.84 

re = T0(-) (1) 

"nth = <th [j^ ) exp(-2: /2p ) , 

where p is the cylindrical radius relative to the 

black hole spin axis, and the assumed constants 
are listed in Table 1 for three fiducial black hole 
spins. In all models, the non-thermal electrons 
have a spectral index of 1.25 and a minimum 
Lorentz factor of 100. The magnetic field strength 
is set to be a fixed fraction (30%) of equipartition 

At frequencies significantly larger than the plasma and cy- 
clotron frequencies, treatments that assume the plasma 
eigenmodes are uncoupled (e.g. Bromley et al. 2001). are 
inappropriate. In the vacuum limit, the two plasma eigen- 
modes remain strongly coupled such that the polarization is 
pajrallely propagated regardless of the details of the plasma. 



relative to the protons and is chosen to be toroidal, 
as found in recent general-relativistic magnetohy- 
drodynamic simulations (De Villiers et al. 2003). 
The accreting gas is assumed to be in free fall in- 
side of the ISCO (6M, 4.233M and 1.237M for 
a = 0, 0.5 and 0.998, respectively), and in Ke- 
plerian rotation otherwise. In all cases the disk 
angular momentum was aligned with the spin of 
the black hole. 

3. Discussion 

As shown in Figure 1, the photosphere be- 
comes smaller at higher frequencies as expected. 
However, special relativistic aberration and the 
Doppler effect make the photosphere asymmetric 
about the rotation axis of the accretion flow. At 
high frequencies, the emission is dominated by a 
small region located on the approaching side of 
the disk. This has two immediate consequences: 
a high polarization at high frequencies (10-20%), 
and a definite shift in the image centroid as a func- 
tion of black hole spin due to strong gravitational 
lensing. At frequencies > 340 GHz the black hole 
"shadow" discussed by Falcke et al. (2000) (see 
also Takahashi 2004; Miyoshi et al. 2004) becomes 
apparent. However, in contrast to previous work, 
the rotation of the disk makes the shadow highly 
asymmetric and even for the non-rotating case it 
only roughly follows the photon capture cross sec- 
tion (shown by the white dashed line). Because 
rotation is a generic feature of the accretion flow, 
it is unlikely that sub-millimeter imaging of Sgr 
A* will show a well defined shadow, as suggested 
before based on a much more simplistic model for 
emissivity profile. 

The unpolarized and polarized spectra are 
shown in Figure 2. High black hole spin (and thus 
high Keplerian velocity at the disk inner edge) 
leads to a broadened sub-millimeter bump, yield- 
ing a marginally better fit to the data. Given 
the current uncertainties inherent in modeling 
the inner edge of the accretion flow, this is not 
necessarily evidence for spin. However, it does 
demonstrate the importance of relativistic effects 
between 30 and 300 GHz. The polarized spec- 
tra shows the anticipated asymptotic behavior at 
high observing frequencies, which strongly dis- 
criminates between low and high spin black holes. 
In addition to the high frequency behavior there is 
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Fig. 1. — Images of an accretion disk around a non- rotating (bottom), moderately- rotating (middle), and 
maximally-rotating (top) black hole at 150 GHz (left), 340 GHz (middle) and 1000 GHz (right) viewed at an 
angle of 45° relative to the spin axis. Polarization tickmarks are overlayed, the lengths of which are indicative 
of the amount of polarized intensity (not polarization fraction). The bar in the upper-right corner indicates 
the scale 2M 10/ias). For reference, the photon capture radius of a non-rotating black hole is delineated 
by the dashed white line. The intensity scale (red being bright and blue being faint) is normalized separately 
for each image and the vertical axis is aligned with the black hole spin. The increasing asymmetry in the 
images is primarily due to the special relativistic effects, which become more significant at the ISCO for a 
rapidly rotating black hole. The dim region near the center of the 150 GHz images is due to the funnel region 
of the thick disk, not due to the black hole "shadow" . The high frequency emission, arising from closest to 
the black hole, is significantly offset between a = ai^ a = 0.998. 



Table 1 

Accretion Model Parameters 



a{M) 


(cm--') 


r," (K) 




Pulh 





3 X 10^ 


1.7 X 10" 


8 X 10* 


-2.9 


0.5 


3 X 10^ 


1.4 X 10" 


5 X 10^ 


-2.8 


0.998 


1 X 10^ 


1.5 X 10" 


1 X 10^ 


-2.8 



'^In the following figures, the short-dash, long-dash and 
soUd lines correspond to the a = 0, a = 0.5 and a = 0.998 
models, respectively. 



a sub-millimeter bump in the polarization fraction 
as well. Again the data appears to marginally fa- 
vor high spin. However, in this case as well, there 
is significant uncertainty in the emission character- 
istics, most notably the geometry of the magnetic 
field. Nevertheless, our assumption of a toroidal 
field is likely to be justified in the inner regions of 
the disk due to the strong shear present in that 
region. 

The position of the image ccntroids as a func- 
tion of frequency is plotted in Figure 3. As sug- 
gested by Figure 1 the image centroids asymp- 
tote to differing fixed positions at high frequen- 
cies. The asymptotic positions are determined by 
strong gravitational lensing, and thus are sensitive 
to the parameters of the black hole. In contrast, 
when the photosphere has a large radius (i.e., at 
low frequencies), the black hole spin has little in- 
fluence and thus the image ccntroids become de- 
generate. The centroid location rapidly changes as 
the photosphere shrinks to the scale of the hori- 
zon. This necessarily occurs during the transi- 
tion between the optically thick and the optically 
thin regimes. Therefore, high precision absolute 
astromctry of the image ccntroids at radio and in- 
frared wavelengths can constrain the black hole 
spin. 

While quantitative differences exist for differ- 
ent magnetic field geometries in both the polarized 
spectra and the locations of the image centroids, 
the qualitative characteristics (which depend only 
upon the asymmetric opacity and relativistic ef- 
fects) remain the same. As a result, detailed multi- 
wavelength studies of the polarization and image 
centroids provide a method by which the spin of 



the black hole at the Galactic center may be mea- 
sured. 
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Fig. 2. — The spectral flux (upper panel) and lin- 
ear polarization fraction (lower panel) are shown 
for the non-rotating {a — 0, short-dash), moder- 
ately rotating (a ~ 0.5, long-dash) and maximally 
rotating (a = 0.998, solid) black hole. The ac- 
cretion disk models are listed in Table 1 and the 
disk is viewed at an angle of 45°. In all cases 
the spectra give an acceptable fit to the existing 
data with the maximally rotating black hole model 
giving the best fit. The polarization fraction ap- 
pears to be a sensitive diagnostic of the spin of 
the black hole near 350 GHz and above 10^ GHz 
(30/Ltm). Below 100 GHz it is believed that Sgr 
A* is Faraday depolarized, and thus our calcula- 
tion (in which Faraday rotation was ignored) is 
inapplicable. The data are taken from Aitken 
et al. (2000); Bower et al. (2002, 2003). The po- 
larization at high frequencies is a direct result of 
the localization of the emission evident in Figure 
1. 
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Fig. 3. — The position of the centroid of the 
disk image is shown as a function of observing 
frequency (labeled in GHz) for the three models 
listed in Table 1 viewed at an angle of 45°. The 
position shifts are measured relative to the black 
hole (where the y-axis is aligned with the black 
hole spin). The projection of the trajectories on 
the Ax axis are shown as functions of frequency 
in the inset. For Sgr A*, the angular scale of the 
axes is M ~ 5/uas times the black hole mass in 
units of 4 X lO^Mfl. 
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